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Abstract 

The experimental density dependence of gas (argon and nitrogen) permeabil- 
ity of partially densified silica aerogels in the Knudsen regime is quantitatively 
accounted for by a computer model. The model simulates both the structure 
of the sintered material and the random ballistic motion of a point particle 
inside its voids. The same model is also able to account for the density de- 
pendence of the specific pore surface as measured from nitrogen adsorption 
experiments. 

PACS numbers: 05.60,+w, 66.30.-h, 61.43.-j, 47.55.Mh 
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The relation between fluid transport properties in porous materials and their geometrical 
characteristics is the subject of great experimental and theoretical activities, expecially after 
the new insights brought by modern scaling ideas and computer modelling facilities |1],|2[]. In 
the case of gas transport in the Knudsen regime, where the molecular diffusion is dominated 
by collisions with the pore walls, it has been demonstrated that the diffusion properties 
strongly depend on the geometrical characteristics of the pore network ||. Therefore, to 
properly account for any experimental results there is a need for a right modelization of the 
internal structure of the material on a mesoscopic scale. 

In this letter we present a comparison between experimental data and numerical cal- 
culations for the gas permeability and the specific pore surface of partially densified silica 
aerogels. The original (non-densified) aerogel structure is modelized by the diffusion-limited 
cluster-cluster aggregation model j|] which has been recently shown to successfully explain 
small angle scattering experiments on these materials 0. The densification process is de- 
scribed by using a numerical realisation of recent scaling ideas for sintering |J and the gas 
diffusivity is calculated by allowing a point particle to move within the pores of the sintered 
structure along successive random straight lines. 

To build the original (non sintered) aerogel structure, we have considered a three dimen- 
sional off-lattice cluster-cluster model J7|!f which has already been extensively described 
elsewhere [f|. A set of N identical spherical particles of diameter a are initially randomly 
disposed (without overlap) in a cubic box of edge length La , occupying a volume fraction, 
or dimension-less concentration of Cq = (ir/6)(N/ L 3 ). In the following the particles are 
assumed to have the density ps (=2.2 g/cm 3 ) of pure silica so that the actual density is 
Cops- Then, these particles are allowed to undergo a brownian diffusive motion and they 
irreversibly stick when they come at contact. Aggregates of particles are also able to diffuse 
together with the individual particles and to stick to particles or to other aggregates. In 
this process, the diffusion constant of the aggregates is considered to vary as the inverse of 
their radius of gyration and periodic boundary conditions are assumed at the box edges. 
When the concentration is sufficiently large (larger than a threshold value c g which tends 
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to zero for infinite box size ||), the final structure is a set of tangent spheres which can be 
described as a loose random packing of connected fractal aggregates, of fractal dimension 
D ~ 2, whose mean size £o decreases as Co increases. 

The sintering process, which consists in a strenghening of the internal structure accom- 
panied by a gradual elimination of the pores, is here numerically simulated by following 
closely the reasoning of reference [f|. Introducing a sintering parameter s, which varies from 
zero (for the non sintered material) to a co-dependent upper limit (for the compact situation 
without hole), we first consider a s- dependent dressed structure in which each initial sphere 
of diameter ao is replaced by a sphere of same center but of larger diameter ad = a (l + s). 
Then the total volume Vd(s) located inside the overlapping spheres is numerically calculated 
(avoiding multiple countings of overlaps). 

The sintered structure is obtained from the dressed structure after applying a proper 
homogeneous volume contraction which insures mass conservation. The length reduction 
factor P(s) is given by: 



is the volume fraction of the sintered structure. 

In figure 1 we provide a typical example with Co = 0.1136. A cut through the initial 
structure is shown as well as at two different stages (s = 1 and s — 2) of the sintering process. 
Although in the last stages (when c(s) becomes close to unity) the shapes of the holes are 
probably not well reproduced (they are certainly less smooth than in reality), this numerical 
method gives a quite good variation of the lower cut-off a as a function of the aerogel density 
p = c(s)ps, as it can be obtained from low angle neutron scattering experiments on partially 
densified silica aerogels. A detailed presentation of this numerical procedure, comparison 
with other numerical methods and scaling predictions and confrontation with experiments 
will be reported elsewhere [|]. 




(1) 



where 




(2) 
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Here, we would like to focus on another geometrical characteristics which is measurable 
by adsorption experiments, namely the specific pore surface area. Knowing that the surface 
of the dressed structure is simply equal to 2dVd/dad = (2/ao)(dVd(s)/ds), the specific surface 
S(s) of the sintered material can be obtained after correcting by the adequate scaling factor 
and dividing by the total mass: 

w n 2 1 dV d 2 dlogc(s) 

a>oPsV d {s) ds a p s ds 

Starting from a given non-sintered structure, characterized by L and cq, one can calculate 
c(s) (through Vd(s)) as well as its derivative, and, using formula (3), one can calculate 
the specific pore surface cir<3cl clS db function of s. Some results for the variation of the 
dimensionless quantity a ps£(s) as a function of c(s) are reported on the log-log plots of 
figure 2 (dashed curves). One has checked that, for s = 0, one recovers the exact result 
aoPs^O) = 6 (horizontal straight line shown in the figure), which is due to the fact that 
in that case the structure is made of identical non overlapping spheres, each of surface area 
7raQ and mass ps^a^. Here the results have been obtained with L = 40.4 and they result 
from an average over 20 realisations of the aerogel structure for each Co value. Here and in 
the following we have checked that the results are almost independent on L. It should be 
noticed that the scaling predicted in ref. f| (a log-log slope of .0/3(3 — D) ~ —0.7) is here 
only observed in a quite small intermediate-c region indicating that corrections to scaling 
are important both at low c (c ~ cq) and at large c (c ~ 1). Obviously £ tends to zero as c 
tends to one. 

To calculate the gas permeability K(s), which is the ratio between a gas flux through the 
sintered material and the applied pressure gradient, we have used the fact that it is simply 
related to the gas diffusion constant D s (s) inside the pores by: 

K(s) = (l-c(s))D s (s) (4) 

Introducing the gas diffusion constant inside the pores of the dressed structure D d (s) = 
j3(s) 2 D s (s) this formula also writes: 



Since the experimental situation that we would like to describe correspond to low pressures 
and small gas particles (smaller than silica particles), the diffusion constant Dd(s) is related 
to the Knudsen motion of a point particle inside the pores of the dressed structure. Given 
the dressed structure, characterized by L, c and s, we have indirectly computed D d (s) by 
allowing a hard sphere of diameter sa to perform a Knudsen motion inside the voids of the 



undressed (non sintered) structure made of spheres of diameter ao- As already noticed |T0 
this is strictly equivalent to the original problem. 

Initially a sphere of diameter sa is released at a random point under the condition 
that it should not overlap any other sphere (when s and/or cq are too large, many trials 
are necessary) and a random direction (uniformly distributed in space) is selected. Then 
the moving sphere is allowed to follow a random straight line motion in this direction until 
it collides a sphere of diameter ao- Immediately after collision a new random direction is 
selected (in a half space) according to the so-called Knudsen cosine law [f^UTl], i. e. such 
that its probability distribution is proportional to cos 9d 2 Q where 9 is the angle between 
this direction and the normal (center-to-center direction) and where d 2 Q = sin OdOdcf) is the 
elementary solid angle. In practice 9 is chosen such that cos 2 9 is uniformly distributed 
between zero and one while the azimutal angle is chosen uniformly between and 2tt. It 
has been shown that this is essential to recover the equilibrium Boltzmann statistics and 
the right transport coefficients Then, after a large number N c of collisions and using 

the periodic boundary conditions, one calculates the end-to-end square displacement £ 2 of 
the diffusing sphere center as well as the total legth A of its trajectory, which is the sum of 
the lengths of the successive segments. This allows us to calculate Dd(s) by the well-known 
formula: 

D d (s) = (6) 



where v = ^j8kT/iim is the mean molecular velocity of the gas particle of mass m at 
temperature T. Reporting in formula (5), one finally gets: 



, 1 1 - c(s) C- 

Some results for the variation with c(s) of the dimensionless quantity K(s)/a v are reported 
in the log-log plot of figure 3. In such calculation we have taken A ~ 10L and we have 
averaged over 10 realizations of the original structure as well as 60 trials for the initial 
position, for each value of c . As expected K tends to zero as c tends to one. The curves 
obtained for different Cq are quite close to each other, showing that the density dependence of 
the permeability is not strongly sensitive to the details of the structure (tangent or overlaping 
spheres). However we have shown above (see figure 2) that considering overlaping spheres 
to simulate the sintering is essential to account for a decrease of the specific pore surface 
with increasing densities. 

Our simulation also allows us to investigate the statistics of the chords, which are the 
segments between successive collisions. It is known that the chord length distribution enters 
many properties of a porous material [j3,l^|. In particular the mean chord length < I > is 



directly related to the specific pore surface via a well known mathematical theorem |T_3[ . In 
our simulation the mean chord length of the dressed structure can be simply obtained by 
< £ >= -if-. Then going back to the sintered structure and using the theorem, we can derive 
the following alternative expression for E(s): 

E( g )=4 1 -f ) (8) 
c(s) p s <£> 

We have reported a ps^ as a function of c(s), calculated from < £ > with the same c values 
as before on the same plot of figure 2. While the overall c-dependence is qualitatively the 
same, we observe some discrepancies compared to the preceeding calculation which increase 
for larger cq. Even if we have checked that the discrepancies do not result from introducing 
an upper cut off to the chord length at A, the numerical uncertainties are certainly larger in 
the second calculation. It may also happen that, for larger concentrations, the presence of 
closed pores may affect differently the two calculations (closed pore surface is fully considered 
in the first calculation). Moreover, the theorem holds for a globally isotropic pore network 



and, even if all the cluster-cluster growth conditions are isotropic, the resulting gel may be 
not since it is known that the cluster-cluster aggregates are intrinsically anisotropic ||14|| . 

The gas permeability experiments on partially densified aerogels have been performed 
using a newly designed apparatus which allows to measures simultaneously the gas flux 
through a cylindrical sample and the gas pressure difference between the two edges. Tech- 
nical details will be given elsewhere fll5| . Here we report on experiments done with a series 



of aerogels which have been prepared under neutral conditions and for which pore surface 
area data have been previously measured by nitrogen adsorption experiments f|. Before 
sintering, their density is known to be equal to po = 0.25 g/cm 3 . Two series of experiments 
have been done with different gas, nitrogen and argon. The experimental results for both 
the ratio K/v and the pore surface area £ are reported in figure 4. As expected, one ob- 
serves that, after dividing by the gas velocity, the data are superimposed. In both cases the 
symbols correspond to the experimental data and the full lines correspond to a fit with our 
computer model using c = Po/p,s = 0.1136, L = 40.4, and a = 70 A . It is remarkable 
that, with only one adjustable parameter (a ), we are able to fit quite well the two quanti- 
ties simultaneously (formula (3) has been used for £). This gives a good confidence in the 
pertinence of our model. Although of the right order of magnitude, the value chosen for a 
is larger than other previous estimates for the mean particle diameter of neutral aerogels 
||16||. Here, we have considered a set of monodisperse spheres while it is well known that the 



size polydispersity of individual particles is very large in neutral aerogels |]I6 |. As already 



noticed for small angle scattering experiments |17j], the large balls dominate any scattering 



process and they also have a greater contribution to the internal pore surface area. This 
could explain why we are obliged to consider here an effective particle diameter larger than 
the simple mean. Also, counting the closed pores (which are consisered in the fit) could also 
increase the effective a . 

In conclusion, using a computer model which simulates both the formation and den- 
sification process of aerogels we have been able to quantitatively account for the density 
dependence of the gas permeability and pore surface area. In the future we will pursue the 



numerical simulations of aerogels sintering, and we will extend the present model, introduc- 
ing polydispersity and analysing more extensively the chord length distribution. We have 
also in project to simulate the brownian diffusion regime in connection with the transport 
properties of fluids inside aerogels. 
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pou, M. Foret and J. Pelous. One of us (A. H.) would like to acknowledge support from 
CONICIT(Venezuela). Laboratoire de Sciences des Materiaux Vitreux is Unite Associee au 
CNRS No. 1119. 
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FIGURES 

FIG. 1. Cut of the simulated aerogel structure at three stages of sintering, s = 0, 1 and 2, where 
c(s) = 0.1136, 0.54, 0.88, repectively. Here and in the other figures the box size is L = 40.4. 

FIG. 2. Log-log plot of aopsT, versus c. Dotted, dashed and dot-dashed curves correspond to 
formula (3) with cq = 0.05,0.1136 and 0.2, respectively. The horizontal line correspond to the 
exact result ao/OsX = 6 for c = cq. Triangles, diamonds and squares correspond to formula (8) with 
Co = 0.05,0.1136 and 0.2, respectively. 

FIG. 3. Log-log plot of K/aov versus c. Triangles, diamonds and squares correspond to 
Co = 0.05,0.1136 and 0.2, respectively. 

FIG. 4. Fits of both K/v and S for neutral aerogels. In the case of K/v open squares and 
diamonds correspond to argon and nitrogen, respectively. The lines are the fits with L = 40.4, 
c = 0.1136 and a = 70A . 
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